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Abstract 

In this article, we analyze the vertices D*D S K, D*DK, DqD s K and 
D s qDK within the framework of the light-cone QCD sum rules approach in 
an unified way. The strong coupling constants Gd*d s k an d Gd*dk are im- 
portant parameters in evaluating the charmonium absorption cross sections 
in searching for the quark-gluon plasmas, our numerical values of the Gp* d s k 
and Gd*dk are compatible with the existing estimations although somewhat 
smaller, the SU(4) symmetry breaking effects are very large, about 60%. For 
the charmed scalar mesons Dq and D s q, we take the point of view that they 
are the conventional cu and cs mesons respectively, and calculate the strong 
coupling constants Gd q d s k and Gd s0 dk with the vector interpolating cur- 
rents. The numerical values of the scalar- D S K and -DK coupling constants 
Gd d 3 k and Gd s0 dk are compatible with the existing estimations, the large 
values support the hadronic dressing mechanism. Furthermore, we study the 
dependence of the four strong coupling constants Gd*d 3 k, Gd*dk, Gd d s k 
and Gd s0 dk on the non-perturbative parameter 04 of the twist- 2 K meson 
light-cone distribution amplitude. 

PACS numbers: 12.38.Lg; 13.25.Jx; 14.40.Cs 
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1 Introduction 

The suppression of the J/ip production in relativistic heavy ion collisions maybe 
one of the important signatures to identify the possible phase transition to the 
quark-gluon plasma [Tj. The dissociation of the J/ip in the quark-gluon plasma 
due to color screening can lead to a reduction of its production, however, the J /if) 
suppression maybe already present in the hadron-nucleus collisions. It is necessary 
to separate the absorption of the J/ip by the nucleons and by the co- mover light 
mesons (ir, K, p, u>, etc.) before we can make a definitive conclusion about the 
formation of the quark-gluon plasma. It is of great importance to understand the 
J/ijj production and absorption mechanisms in the hadronic matter. The values of 
the J/ip absorption cross sections by the light hadrons are not known empirically, 
we have to resort to some theoretical approaches. Among existing approaches for 
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evaluating the charmonium absorption cross sections by the light hadrons, the one- 
meson exchange model and the effective £77(4) theory are typical [21 El- The detailed 
knowledge about the hadronic vertices or the strong coupling constants which are 
basic parameters in the effective Lagrangians is of great importance. 

The discovery of the two strange-charmed mesons D s o and D s \ with spin-parity 
+ and 1 + respectively has triggered hot debate on their nature, under-structures 
and whether it is necessary to introduce the exotic states jl]. The mass of the D s0 
is significantly lower than the values of the + state mass from the quark models 
and lattice simulations [5 . The difficulties to identify the D s q and D s \ states with 
the conventional cs mesons are rather similar to those appearing in the light scalar 
mesons below lGeV. Among the various explanations, the hadronic dressing mech- 
anism is typical. The scalar mesons a (980), / (980), D and D s0 may have bare qq, 
cu and cs kernels in the P— wave states with strong coupling to the nearby threshold 
respectively, the S— wave virtual intermediate hadronic states (or the virtual mesons 
loops) play a crucial role in the composition of those bound states (or resonances 
due to the masses below or above the thresholds). The hadronic dressing mechanism 
(or unitarized quark models) takes the point of view that the /o(980), ao(980) , D 
and D s o mesons have small qq, cu and cs kernels of the typical qq, cu and cs mesons 
size respectively. The strong couplings to the virtual intermediate hadronic states 
(or the virtual mesons loops) may result in smaller masses than the conventional 
scalar qq , cu and cs mesons in the constituent quark models, enrich the pure qq , cu 
and cs states with other components (HI Ej . Those mesons may spend part (or most 
part) of their lifetime as virtual KK, D S K and DK states (Eld- It is interesting to 
study the possibility of the hadronic dressing mechanism. 

In this article, we calculate the values of the strong coupling constants Gd*d 3 k, 
Gd*dk, Gd d s k and Gd s0 dk within the framework of the light-cone QCD sum rules 
approach. The light-cone QCD sum rules approach carries out the operator prod- 
uct expansion near the light-cone x 2 ~ instead of the short distance x ~ while 
the non-perturbative matrix elements are parameterized by the light-cone distri- 
bution amplitudes which classified according to their twists instead of the vacuum 
condensates [SI El- Furthermore, we study the dependence of the strong coupling 
constants Gd*d s k, Gd*dk, Gd d s k and Gd s0 dk on the coefficient 04 of the twist- 2 
K meson light-cone distribution amplitude 4>k{u), and estimate the values of the 
non-perturbative parameter. It is very difficult to determine the with the QCD 
sum rules, the values of the 04 suffer from large uncertainties, as it concerns high 
dimension vacuum condensates which are known poorly [HJ El El E] • It is of 
great importance to determine the values directly from the experimental data. 

The article is arranged as: in Section 2, we derive the strong coupling constants 
Gd*d s k, Gd*dk, Gd d s k and Gd s0 dk within the framework of the light-cone QCD 
sum rules approach; in Section 3, the numerical results and discussions; and in 
Section 4, conclusion. 
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2 Strong coupling constants G D * DsK , G D * DK , G DqDsK 
and Gn s0 DK with light-cone QCD sum rules 

In the following, we write down the definitions for the strong coupling constants 
Gd*d s k, Gd*dk, Gd d s k and Gn a oDK i 

(D*(q + P)D s (q)\K(P)) = G D « DaK (P - q) ■ e , 

(D* s (q + P)D{q)\K(P)) = G D . DK (P ~ q) • e , 

(D (q + P)D s (q)\K(P)) = G DoDaK , 

(D s0 (q + P)D(q)\K(P)) = G DaoDK , (1) 

here the e M are the polarization vectors of the mesons D* and D*. We study the 
strong coupling constants Gd*d s k, Gd*dk, Go D a K and Gd s0 dk with the interpo- 
lating currents Jd s (x), Jd(%), J? s (x) and J^{x) in an unified way, and choose the 
two-point correlation functions I1^(P, g) and IT 2 (P, q), 
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= c(x)i7 5 s(x) . 






(4) 



The correlation functions 11^ (P, g) can be decomposed as 

nj (2) (^) = 4 (2) (g 2 ,(g + P) 2 )P M + nf)(g 2 ,(g + P) 2 )g„ (5) 

due to the Lorentz covariance. In this article, we derive the sum rules with the 
tensor structures P M and g M respectively, and make detailed studies. 

According to the basic assumption of current-hadron duality in the QCD sum 
rules approach ^3], we can insert a complete series of intermediate states with the 
same quantum numbers as the current operators Jo 3 (a;) (Jd(x)) and J?(x) (Ju s ( x )) 
into the correlation function IT* (II 2 ) to obtain the hadronic representation. After 
isolating the ground states and the first orbital excited states contributions from 
the pole terms of the D s , D* and D (D, D* and D s0 ) mesons, the correlation 
function IT* (II 2 ) can be expressed in terms of the strong coupling constants G and 
the decay constants /a/ of the heavy mesons, the explicit expressions are presented 
in the appendix. We use the standard definitions for the decay constants fu (/r» s) 
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Id, /dj, Id*, /d. , /a> ) of the heavy mesons, 



<0|Jj,(0)|D(g)> 



(0|J D .(0)|D a (g)> 



fD a m 2 Da 



(0\J°(0)\D*(q)) 
(0|J M D (0)|A)(<Z)> 

(o|Jf(o)|^o(g)) 



f D: m D; e li , 

fDoQ/j. j 




(6) 



The quarks c and s have finite and non-equal masses, the non-conservation of 
the vector currents J^ s (x) and J®(x) can lead to the non- vanishing couplings to the 
scalar mesons D s q and Dq beside the vector mesons D* and D*, we can study the 
properties of those mesons with the two interpolating currents J® s (x) and JF(x) in 
an unified way. Here we have not shown the contributions from the high resonances 
and continuum states explicitly as they are suppressed due to the double Borel 
transformation. The numerical values of the fractions 



are less than 30% and the corresponding spectral densities for the ground states are 
greatly suppressed, the tensor structures with are especially suitable for studying 
the first orbital excited states D and D s0 with the vector currents. The numerical 
values of the fractions 



are about 2, the tensor structures with P M are especially suitable for studying the 
ground states D* and D* with the vector currents. 

Now we carry out the operator product expansion near the light-cone x 2 ~ 
to obtain the representation at the level of quark-gluon degrees of freedom for the 
correlation functions and IT^. In the following, we briefly outline the operator 
product expansion for the correlation functions 11^ and in perturbative QCD 
theory. The calculations are performed at the large space-like momentum regions 
(q + P) 2 <^ and q 2 <^ 0, which correspond to the small light-cone distance x 2 ~ 
required by the validity of the operator product expansion approach. We write down 
the propagator of a massive quark in the external gluon field in the Fock-Schwinger 




2 i 2 ™2 -_,2 _i_ _2 

m D a + m K m D* s -m D + m K 

2 ' 2 



m 2 D » + m 2 Ds - m 2 K m 2 D * + m 2 D - m\ 



2 ' 2 
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gauge firstly [TUj . 



(o|r{ ft (x 1 )^)}|o> = i J 

1 

kP-m 2 ^ 3 ~ I dV9s G %( VXl + V ^ X2) 
o 

1 Jt + m 1 l\ 

here is the gluonic field strength, g s denotes the strong coupling constant. 
Substituting the above c quark propagator and the corresponding K meson light- 
cone distribution amplitudes into the correlation functions IT* and 11^ in Eqs.(2- 
3) and completing the integrals over the variables x and k, finally we obtain the 
representation at the level of quark-gluon degrees of freedom, the explicit expressions 
are presented in the appendix. In calculation, we have used the two-particle and 
three-particle K meson light-cone distribution amplitudes [HI El EEH HU El- the 
explicit expressions are also presented in the appendix. The parameters in the light- 
cone distribution amplitudes are scale dependent and can be estimated with the 
QCD sum rules approach [SJ El US EEU U2]- In this article, the energy scale /i is 
chosen to be /i = lGeV. 

We perform the double Borel transformation with respect to the variables Q\ = 
— (q + P) 2 and Q\ = —q 2 for the correlation functions Hp and H^ 2 \ and obtain 
the analytical expressions for those invariant functions, the explicit expressions are 
presented in the appendix. 

In order to match the duality regions below the thresholds So and s' for the 
interpolating currents J^(x)(J^ s (x)) and Jd s (x) (Jd( x ) ) respectively, we can ex- 

1(2) 1(2) 

press the correlation functions n p and IV at the level of quark-gluon degrees of 
freedom into the following form, 



P(s, s') 
{s - {q + P) 2 } {s> - q 2 ) 



4 ( g)(g 2 ,(g + P) 2 ) = Ids Ids' - ' (8) 



then we perform the double Borel transformation with respect to the variables Q\ = 
— (q+P) 2 and Q2 = —q 2 directly. However, the analytical expressions for the spectral 
densities p(s, s') are hard to obtain, we have to resort to some approximations. As 
the contributions from the higher twist terms are suppressed by more powers of -^p 
or _ ( - g ^p- )2 ; the continuum subtractions will not affect the results remarkably, here 
we will use the expressions in Eqs. (28-29) for the three-particle (quark-antiquark- 
gluon) twist-3, twist-4 terms, and the two-particle twist-4 terms. In fact, their 
contributions are of minor importance, the dominating contributions come from the 
two-particle twist-2 and twist-3 terms involving the <Pk{u), 4> p (u) and (T ('u). We 
perform the same trick as Refs. (TUl EI and expand the amplitudes 4>k(u), 4> p {u) and 
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(u) in terms of polynomials of 1 — u 

) d 

' du 



N 



k 



2 

m 



k=0 
N 

then introduce the variable s' and the spectral densities are obtained. After straight- 
forward but cumbersome calculations, we can obtain the final expressions for the 

1(2) 

double Borel transformed correlation functions n M at the level of quark-gluon 
degrees of freedom below the thresholds. The masses of the charmed mesons are 
Mb* = 2.012GeV, M D . = 2.112GeV, M D = 1.865G?eV, M Ds = l.97GeV, M Do = 
2A0GeV, and M DsQ = 2.317GW, the ratios are Md +^ d - « 0.45, AfD +g fg . ~ 0.47, 

Mr^+M D * ~ 0.49 and ^P+aS ~ 0.45 [T3]. There exist overlapping working win- 
dows for the two Borel parameters M\ and M|. It's convenient to take the value 
M\ = Ml u = j^fe, = |, M 2 = = \Ml furthermore, the K meson 

light-cone distribution amplitudes are known quite well at the value it = § compar- 
ing with the values at the end-points. We can introduce the threshold parameter so 
and make the simple replacement, 

f m\ + tt (l - u )m 2 K \ f m\ + u (l - u )m 2 K \ f s } 

eXP \ M 2 j^ ex Pj M2 | exp| M2 | 

to subtract the contributions from the higher resonances and continuum states ^Hj , 
finally we obtain the following sum rules, 

G D * DsK m D ,f D ,f Da m 2 D m 2 D ,+m 2 D -m 2 K f m 2 D , m 2 D \ 
m~^ s ml, eXP l-Mf- -Mj\ = AA > 



G p*pK m D* fp* fp m p m Dt +m D -m z K J m D% m D 
m c + m u m 2 D , 6XP | Ml M? 



G D *D s Krn D *f D *f Da m 2 D3 m D *-m Dg +m 2 K ( m D * m 



exp 



m c + m s m 2 D , y Ml M; 

+ G DoDsK f D J D m 2 Ds ^ _ < | = 



m c + m s { Ml M. 



G DtDK m D *f D *f D m 2 D m 2 Dt - m 2 D + m 2 K j m 2 D% m \ 
m c + m u m 2 D * e P j ^2 ^2 

G DsoDK f D J D m 2 D |_m^ _ (13) 
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G D * DsK m D *f D *f Ds m 2 Ds m 2 D , + m 2 Ds 
m c + m s m 2 D * 
G d d s k fD /d s m 1 s 



exp 




(14) 



m c + m, 



Go* DKm Dt f Dt f D m 2 D m 2 Dt + m 2 D - m\ 




dkIdmId^d 





(15) 



The explicit expressions of the notations AA, BB, CC, DD, EE and FF are 
lengthy and given explicitly in the appendix. A slight different manipulation (with 
the techniques taken in the Ref. [T%l IT!?] ) for the dominating contributions come 
from the terms involving the two-particle twist- 2 and twist-3 light-cone distribution 
amplitudes 4>k(u), 4> p (u) and 4> a (u) leads to the sum rules with the same type as in 
Ref . [IH] • However, those type sum rules are not stable with respect to the variations 
of the Borel parameter M 2 , here we will not show the expressions explicitly for 
simplicity. It is not surprise that the QCD sum rules as a QCD model have both 
advantages and shortcomings. 

3 Numerical results and discussions 

The input parameters are taken as m s = (140 ± 10)MeV, m c = (1.25 ± 0A0)GeV, 
A 3 = 1.6 ± 0.4, f 3K = (0.45 ± 0.15) x 10" 2 GW 2 , u 3 = -1.2 ± 0.7, u 4 = 0.2 ± 0.1, 
a 2 = 0.25±0.15, ai = 0.06±0.03 [H] 13 GDI EED EES , fx = 0A60GeV, m K = A98MeV, 
m Ds0 = 2.317GeV, m D = 1.865GeV, = 1.97GeV. In this article, we take the 
values of the a 4 to be zero, and explore the dependence of the strong coupling 
constants Gd*d s k, Gd*dk, Gn D a K and Gd s0 dk on this parameter. 

For the threshold parameter s° Do , we can use the experimental data as a guide, 
m Do = 2A0GeV, T mDQ = 283MeV [J3], and choose the values s° Dq = (Q.8-7.2)GeV 2 
to subtract the contributions from the high resonances and continuum states. The 
mass and width of the D from Belle and Focus are vtidq — 2308 ± 17 ± 15 ±28MeV, 
T Do = 276 ± 21 ± 18 ± 60MeV [JB], m Do = 2407 ± 21 ± 35MeV, T Do = 240 ± 55 ± 
59MeV ^7j. The predictions from the constituent quark models are mo = 2AGeV 
The values of the mass from the two collaborations have the difference about 
lOOMeV, in this article, we take the value m£> = 2AGeV as input parameter, 
our final numerical results for the large strong coupling constant Gd d s k support 
smaller values for the D if the same mechanism takes place for both the charmed 
scalar mesons D and D s0 . Furthermore, the strong coupling constant Gd d s k 
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* s°=4.7GeV 2 ; 

* s°=4.9GeV 2 ; 

* s°=5.1GeV 2 . 




3.0 3.5 4.0 4.5 5.0 5.5 

M 2 (GeV 2 ) 



Figure 1: The Gd*d s k with the parameters M 2 and s° D , from Eq.(10). 



is not sensitive to the values of the mo , taking the values mo = 2AGeV or 
m Do = 2.3GeV can not change the conclusion qualitatively or quantitatively. 

For the threshold parameters s° D * , s° D * and s° Dg0 , the experimental values of the 
masses are wid* = 2.QlGeV , mo* = 2.\\2GeV and mD a0 = 2.317GeV , the widths 
are very narrow jT^j. We can choose the values of the threshold parameters s° D * = 
(4.7 - 5.1)GeV 2 , s D . = (4.8 - 5.2)GeV 2 and s° DsQ = (7.0 - 7A)GeV 2 to subtract 
the contributions from the high resonances and continuum states. From Figs. 1-3, 
we can see that the numerical values of the strong coupling constants Gd*d s k and 
Gd d b k are not sensitive to the threshold parameters s° in those regions, the values 
we chosen here are reasonable. 

The values of the decay constants fn, /d s) /d*, Id*, /d and fr, s0 vary in a 
large range, for example, f D = (0.17 ± 0.01)Ge^, f D * '= (0.24 ± 0.02)GW [TDJ, 
f Do = (0.217 ±0.025)GeK, m Do = 2.212GeV [20], fD sQ = (0.225 ± 0.025)Ge1/ EH, 
f D = (0.177 ±0.021)GeV, f Ds = (0.205 ± 0.022)GeV [221, f Do = (0.17 ± 0M)GeV 
[231 from the Q CD sum rules ; fD a = 0.268GeV, f D , = 0.315GeV, f D = 0.2MGeV, 
f D . = OMOGeV |21, f D * = 0.375 ± 0.02AGeV, °f D * = 0.340 ± 0.023Ge1/ j2S], 
f D = 0.238Ge1/, f Ds = OMlGeV [26 from the potential models; f D , = 32^3,MeV 



-17 J 



f D . = 223tHMeV j2Z] from the quark models, and f D = (222.6 ± 16.7?£j)MeV 



from the experimental data [28]. For a review of the values of the decay constants 
for the mesons D and D s from the QCD sum rules and lattice QCD, one can consult 
the second article of the Ref. p. 

In this article, we take the following constraints for the decay constants, 



1.0 < 



^f_ 

fo f 



fDs 

fD 



< 1.1 



(16) 
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Figure 2: The Gd d s k with the parameters M 2 and s^ from Eq.(12). 
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Figure 3: The G£> d s k with the parameters M 2 and s^ from Eq.(14). 
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and choose the values, 

f Ds = (0.25 ± 0.02)GeV , f D = (0.23 ± 0.02)GeV, 
f Ds0 = (0.225 ± 0.025)Ge\/ , f Do = (0.217 ± 0.020)GeV , 

/d ; = (0.26 ± 0.02)GeV , f D . = (0.24 ± 0.02)GeV . (17) 

In numerical calculation, we observe that the values of the strong coupling constants 
Gd*d 3 k, Gd*dk, Gd d s k and Gd s0 dk are sensitive to the six hadronic parameters, 
small variations of those parameters can lead to relatively large changes for the 
numerical values, refining the six hadronic parameters is of great importance. 

The Borel parameters in Eqs.(lO-ll) are taken as M 2 = M| = (6 - 12)GeV 2 
and M 2 = (3 — Q)GeV 2 , in those regions, the values of the strong coupling constants 
Gd*d s k and Gd*dk are rather stable from the sum rules in Eqs. (10-11) with the 
simple subtraction, which are shown, for example, in the Fig.l and Figs. 4-7 for the 
strong coupling constant Gd*d s k, similar figures can be obtained if the values of 
the strong coupling constant Gd*dk are plotted. In this article, we only show the 
numerical values from the sum rules in Eq.(10), Eq.(12) and Eq.(14) explicitly for 
simplicity. 

The Borel parameters in Eqs. (12-15) are chosen as M\ = Mf = (10 — 20)GeV 2 
and M 2 = (5 — 10)GeV 2 , in those regions, the values of the strong coupling con- 
stants Gi) a D s K and Gd s0 dk are rather stable from the sum rules in Eqs. (12-13) with 
the simple subtraction, which are shown in the Fig2, Figs. 4-6 and Fig. 8 for an illus- 
tration. However, the strong coupling constants Gd d s k and Gd s0 dk from the sum 
rules in Eqs. (14-15) have a negative sign comparing with the corresponding ones 
from the sum rules in Eqs. (12-13), and much smaller absolute values. The fractions 

m 2 D * + m 2 Da - m 2 K m 2 D , + m 2 D - m 2 K 

2 ' 2 

are about 2. In the sum rules in Eqs. (10-11), the ground state saturate condition can 
be safely satisfied below the threshold s° D * (s%*). The vector interpolating current 
J?(x) (J^ s (x)) has both non- vanishing couplings to the vector state D* (£)*) and to 
the scalar state Dq (D s q), there are two hadronic states, the ground state D* (D*) 
and the first orbital excited state Dq (D s q) in the channel cu (cs) below the threshold 
s° Dq (s° Ds0 ), the ground states D* and D* are not suppressed due to the factor 2, the 
sum rules in Eqs. (14-15) are not suitable for studying the strong coupling constants 
Gd d s k and Gd s0 dk, our final numerical values support this assumption. We show 
this fact in the Fig. 3 for an illustration. 

We determine the values of the strong coupling constants Gd*d s k and Gd*dk 
from the Eq.(10) and Eq.(ll) respectively, then use those values as the input param- 
eters, and calculate the values of the strong coupling constants Gd d s k and Gd s0 dk 
from the Eqs. (12-15) respectively. 

The uncertainties of the five parameters 004, 0U3, A3, m c and a± can not lead to 
large uncertainties for the numerical values. The main uncertainties come from the 
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• f =0.60*1 0" 2 GeV 2 ; 
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Figure 4: The Gd*d s k and Gd d s k with the parameters M 2 and / 3 x from Eq.(10) 
and Eq.(12) respectively. 

ten parameters f 3K , m s , a 2 , r? 4 , / D , / D ., f Do , f Ds and / Ds0 , small variations 
of those parameters can lead to relatively large changes for the numerical values, 
which are shown in the Figs. 4-8 for an illustration. 

Taking into account all the uncertainties, finally we obtain the numerical results 
for the strong coupling constants, 

Gd*d s k = 2.02lo|g , G DoDs k = 6.5t\'*GeV , 
Gd*dk = l-84i^63 , Gd s0 dk = , (18) 

which are shown in the Figs. 9- 10 respectively. 

The strong coupling constants Gd*d s k, Gd*dk, Gd d s k and Gd 30 dk can be 
related to the parameters g and h in the heavy-light Chiral perturbation theory 

Gsptt = V^smp— — —-j-, 

2^/m P m v 
^VP-k — 7 9 j 

In 

here the S are the heavy scalar mesons with + , the P are the heavy pseudoscalar 
mesons with 0~, the V are the heavy vector mesons with 1~, and the 7r stand for 
the light pseudoscalar mesons. 

The parameter g has been calculated with the light-cone QCD sum rules jH21 EH 
Elj , the quark models [SSI IHEj and extracted from the experimental data |H3 EH] • 
The values vary in a large range, the corresponding values of the strong coupling 
constants Gd*d s k and Gd*dk in the SU(3) limit for the light pseudoscalar mesons 
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-*-f D =0.25GeV; 
• f D *=0.27GeV; 
^ f D *=o.23GeV. 
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Figure 5: The Gd*d s k and Gd d s k with the parameters M 2 and fo s from Eq.(10) 
and Eq.(12) respectively. 
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Figure 6: The Gd*d s k and Gd d s k with the parameters M 2 and 02 from Eq.(10) 
and Eq.(12) respectively. 
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-*-f D .=0.24GeV; 

• f D .=0.26GeV; 

* f D .=0.22GeV. 
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Figure 7: The Gd*d s k with the parameters M 2 and fu* from Eq.(10). 

are listed in the Table. 1. From the table, we can see that our numerical results are 
compatible with the existing estimations, although somewhat smaller. 

The values of the strong coupling constants Gd*d s k and Gd*dk are sensitive to 
the non-perturbative parameter 04, if we take a larger value rather than zero, larger 
values of the Gd*d s k and Gd*dk are obtained. The Gd*d s k and Gd*dk are more 
sensitive to the a 4 comparing with the Gd d s k and Gd s0 dk, which are shown in 
the Fig. 11. In fact, the largest uncertainties come from the uncertainties of the 04, 
they are ideal channels to determine this parameter directly from the experimental 
data. Once the experimental data for the values of the strong coupling constants 
Gd*d s k and Gd*dk are available, powerful constraints can be put on the range of 
the parameter 04. If we take the values from the QCD sum rules as input parameters 
pTTj . Gn*£) a K = 3.02 ±0.14 and Gd%dk = 2.84 ±0.31, very large values of the a 4 are 
obtained. 

The parameter h has been estimated with the light-cone QCD sum rules [3*9*] . 
the quark models [""""J], Adler-Weisberger type sum rules "4~0~] . and extracted from 
the experimental data [H], the values are listed in the Table. 2, from those values 
we can estimate the values of the corresponding strong coupling constants Gd d s k 
and Gd s0 dk m the SU(3) limit for the light pseudoscalar mesons. The value of 
the dimensionless effective coupling constant Y/k = 0.46(9) from Lattice QCD (321 
is somewhat smaller than the values extracted from the experimental data Y/k = 
0.73^24, here the Y is the decay width and the k is the decay momentum. Our 
numerical values Gd d s k = Q-5t.\'^GeV and Gd s0 dk = 5-9-\'lGeV are compatible 
with the existing estimations in Refs. |36| 13*9*1 Ell EL although somewhat smaller 
comparing with the values obtained in Ref . |T*""| with the scalar interpolating current 
for the D s0 meson, and about 2 — 3 times as large as the energy scale Md s0 = 
2.317GeV , and favor the hadronic dressing mechanism. For a short discussion about 
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Figure 8: The Gd d s k with the parameters M 2 and fc , m s , 7/ 4 from Eq.(12). 
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Figure 9: The Gd*d s k(&) and Gd*dk(Jo) with the parameter M 2 from Eq.(lO) and 
Eq.(ll) respectively. 
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Figure 10: The Gd d s k(&) and Gd 30 dk(Jo) with the parameter M 2 from Eq.(12) and 
Eq.(13) respectively. 
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\g\ 


Gd*d s k 


Gd*dk 


Reference 


0.38 ±0.08 


9.5 ±2.0 


9.4 ±2.0 




29 




6.04 ±0.28 


5.68 ±0.62 


IHU 


0.34 ± 0.10 


8.5 ± 2.5 


8.4 ± 2.5 




32 


0.28 


7.0 


6.9 


33 


0.35 ±0.10 


8.7 ±2.5 


8.7 ±2.5 


[34 


0.50 ±0.02 


12.4 ±0.5 


12.4 ±0.5 


[35 


0.61 


15.2 


15.1 


36 


0.59 ±0.07 


14.7 ± 1.7 


14.6 ± 1.7 


[37 


n 07+u.ut) 
u - z ' -0.03 


6 7 +i - 5 

°- ' -0.7 




38 


u - iu -0.05 


4.041};^ 


3.68ii- 


This work 



Table 1: Numerical values of the parameter g, and the corresponding values of the 
strong coupling constants Gd*d 3 k and Gd*dk hi the SU(3) limit. Here we have 
double the values of our numerical results and the ones from Ref.jHI] due to the 
difference between the definitions for the strong coupling constants. 

the hadronic dressing mechanism, one can consult Ref.|19j. or one can consult the 
original literatures for the details jHl U\ ■ 

The large values of the strong coupling constants Gd d s k and Gd s0 dk obvi- 
ously support the hadronic dressing mechanism, the D and D s0 (just like the 
scalar mesons /o(980) and ao(980), see Ref.[TH]) can be taken as having small scalar 
cu and cs kernels of typical meson size with large virtual S-wave D S K and DK 
cloud respectively. In Ref.[30j, the authors analyze the unitarized two- meson scat- 
tering amplitudes from the heavy-light Chiral Lagrangian, and observe that the 
scalar meson D s q appears as the bound state pole with the strong coupling constant 
Gd s0 dk = 10.203CreV. Our numerical results Gd s0 dk = 5.9^\'lGeV are smaller, the 
values of our previously work Gd s0 dk = 9.3^2jGel / with the scalar interpolating 
current are more satisfactory [T9~] . 

4 Conclusions 

In this article, we analyze the vertices D*D S K, D*DK, D D S K and D s0 DK within 
the framework of the light-cone QCD sum rules approach in an unified way. The 
strong coupling constants Gd*d s k and Gd*dk are important parameters in evalu- 
ating the charmonium absorption cross sections in searching for the quark-gluon 
plasmas, our numerical values of the Gd*d s k and Gd*dk are compatible with 
the existing estimations although somewhat smaller, the SU (4) symmetry break- 
ing effects are very large, about 60%, the approximation of the SU{4) symmetry 
Gd*d s k = Gd*dk = 5.0 is not suitable r3|. For the scalar mesons D and D s0 , 
we take the point of view that they are the conventional cu and cs meson respec- 
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\h\ 


G DoDsK {GeV) 


G Ds0 D K {GeV) 


Reference 


0.881-- 




".3-t; 


[19J 






10.203 


[30j 


0.536 


5.7 


5.68 


[36j 


0.52 ±0.17 


5.5 ±1.8 


5.5 ±1.8 


[39j 


< 0.93 


< 9.9 


< 9.86 


[40j 


0.57-0.74 


6.1-7.9 


6.0-7.8 


[41j 


0.611SJJ (or 0.561SS) 




5 9+ LV 


This work 



Table 2: Numerical values of the parameter h, and the corresponding values of the 
strong coupling constants Gd d s k and Gd s0 dk m the SU(3) limit. 
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Figure 11: The G D * DsK (a), G D * DK (b), G DoDsK (c), G DaoDK (d) with the parameters 
M 2 and a 4 from Eq.(10), Eq.(ll), Eq.(12), Eq.(13) respectively . 
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tively, and calculate the strong coupling constants Gd d s k and Gd s0 £,k within the 
framework of the light-cone QCD sum rules approach. The numerical values of 
the scelax-D S K and -DK coupling constants Gd d 3 k and Gd s0 dk are compatible 
with the existing estimations although somewhat smaller, the large values support 
the hadronic dressing mechanism. Just like the scalar mesons /o(980) and ao(980), 
the scalar mesons D and D s0 may have small cu and cs kernels of typical cu and 
cs mesons size respectively. The strong coupling to virtual intermediate hadronic 
states (or the virtual mesons loops) can result in smaller mass than the conventional 
scalar mesons cu and cs in the constituent quark models, enrich the pure states cu 
and cs with other components. The Dq and D s q may spend part (or most part) 
of they lifetime as virtual D S K and DK states. Furthermore, we study the depen- 
dence of the strong coupling constants Gd*d s k and Gd*dk on the non-perturbative 
parameter a 4 of the twist-2 K meson light-cone distribution amplitude. The values 
of the strong coupling constants Gd*d s k and Gd*dk are more sensitive to the 
comparing with the Gd d s k and Gd s0 dk- The largest uncertainties come from the 
uncertainties of the 04, they are the ideal channels to determine the parameter di- 
rectly from the experimental data. Once the experimental data for the values of the 
strong coupling constants Gd*d 3 k and Gd*dk are available, powerful constraints 
can be put on the range of the parameter 04. 



18 



Appendix 

The explicit expressions of the correlation functions II* and nj; in the hadronic 
representation, 



K 



<0\JP\ D*(q + P) >< D*Ds I K >< D s (q)\J Ds | > 
{m 2 D , - (q + P) 2 } (ml a - ?) 
<0 | J» | D (g + P) xDpDs \K><D s (q)\J Ds | 0> 
+ {m 2 Do - (q + P) 2 } (mi s - qi) + 

G D * DsK m D *f D *f Ds m 2 Ds (P - q) ■ ee M 
(m c + m s ) {m 2 D , - (q + P) 2 } (m 2 Ds - q 2 ) 
G Do D s KfD fD s m 2 Ds (q + P)^ 



(m c + m s ) {m 2 Do - (q + P) 2 } (m 2 Ds - q 2 ) 



+ 



G D *D s Km D *f D *f Da m 2 D m 2 D * - m 2 D + m 2 



K 



(m c + m s ) {m 2 D , - (g + P) 2 } (m 2 Ds - q 2 ) m 2 D * 

GD D s KfD fD s m 2 D 



q» + 



+ 



(m c + m s ) {m 2 Do - (q + P) 2 } (m 2 Dg - q 2 ) 

G D *D sK m D *f D *f Ds 7n 2 Ds m? D . + m 2 Ds - m\ 

(m c + m s ) {m 2 D , - (q + P) 2 } (m 2 Ds - q 2 ) m 2 D * 

G D D S K f D f D a m \) a 



(m c + m s ) {m 2 Do - (q + P) 2 } (m 2 Ds - q 2 ) 



P„ + ---, (19) 
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n2 < I tf- \ D* s (q + P) >< D* S D | K X D(q)\J D \ > 



{m^-fo + P) 2 }^-? 2 ) 
<0\J^\ Dsojq + P) >< DsoD I g >< g(g)| Jp | > 

G D * DK m D *J D *f D m D (P - g) ■ ee M 
(m c + m u ) |m|, ; - (q + P) 2 j (m|, - q 2 ) 
G d s0 dk f D s0 f Dm 2 D (q + P) M 



+ 



+ m u ) {m 2 DsQ - (q + P) 2 } (m 2 D - q 2 ) 



+ 



G D * DK mp* f D * f D m 2 D m 2 Di - m 2 D + m 2 - 



K 



(m c + m u ) jm 2 ^ - (q + P) 2 | (m 2 D - q 2 ) 



rn 



_! Gp a0 DKfD a0 fD'm 2 D 

(m c + m u ) {m 2 DaQ -(q + P) 2 } {m 2 D - q 2 ^ ' 



G D * DK m D * f D * fpm 2 D m 2 D * + m 2 D - m\ 

{m c + m u )\m 2 D% - (q + P) 2 } (m 2 , - q 2 ) m Dt 
, G Ds0 p K fp s0 f D m 2 D I 

The explicit expressions of the correlation functions n* and II 2 at the level of 
quark-gluon degrees of freedom, 



LI,' = q„ \ I du 



I' 



m s Jo (q + uP) 2 - m 2 c 

Bit) 



pi l>U 

+m c f K m 2 K du dt— 

Jo Jo {{q + uP) 2 - m l c ) 

+fm m2 K I dv [ da g f da 
Jo Jo Jo 



{[ q + ((l- v )a g + a s )P] 2 -m 2 c y 
$(1 -a- (3, (3, a) 



r-1 r-1 j-a g rl-p 

—Am c f k^Ac I dvv l da„ I df3 I da 

Jo Jo Jo Jo {[ g + (i_ m9 )p] 2 _ m 2}° 

+4m c fK r m A K [ dv [ da g f da s (da- 
Jo Jo Jo Jo 



{[ q +((l-v)a a + a a )P] 2 -miy 
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+ Pa 



fi<m 2 K f L u<j) p (u) 



pi PU 

+ m cfK m2 K / du I dt 
Jo Jo 



m s jo 

,2 /•! 



(q + uP) 2 — m\ 



uB(t) 



{(q + uP) 2 — m 2 Y 



+4^ / #.(«) 



6m 



s JO 
1 



+m c f K / du 



l-u 



>K 



d 

du 



2ml 



(q + uP) 2 - m 2 c [(q + U P) 2 - m 2 ] 



m K m c 



A{u) 



o {(q + uP) 2 -m 2 c 2 [( q + U P) 2 - m 2 ] 



i pi 



+f3Rm K dv da 



o Jo 



l-a„ 



da 



[(1 - v)a g + a s } (2v - 3) T(a u , a g: a s ) 
{[ q + ((l- v )a g + a s )P] 2 -m 2 c } 2 



-2f 3K J dw J da g J r/r..;/-(r., i .r. 7 .r 1> ) — — — : - - 



-AmJ K m K I dvv 

'o Jo 



1 r a a r l -P 
da g I d(3 
Jo Jo 



da 



+Am c fKm] < I dv I da g I da s / da 

Jo Jo Jo Jo + ((1 - v )a g + a s )P] 

— m c f ^m 2 ^ I dv I da g [ da s ^— ^ — — — — 

Jo Jo Jo {[q+ ((l- v )a g + a s )P] 2 - m 2 c ) 



du (q + uP) 2 — m 2 
(1 - va g )<&{\ - a - (3, (3, a) 
{{q + (l-va g )P] 2 -m 2 c } 3 
((1 — v)a g + — a — a g , a g , a) 

T 

(21) 



2 2 

— mi 



2 [ ' 



n M = n ^ ( u i — ► l~u;a u < — ► a s ) . 



(22) 
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The light-cone distribution amplitudes of the K meson, 



< 0\u(0) W s(x)\K(P) > = if K P» due 



< 0\u(0)^g s G a p(vx)s(x)\K(P) > 



2 2 

-iuP-x J ,„ /„.\ , m K X 



<Pk{u) + 



16 



A(u) 



<0\u(0)i~f 5 s(x)\K(P) > = 

<0\u(0)a lxul5 s(x)\K(P)> = 
< 0\u(0)a a ^ 5 g s G^(vx)s(x)\K(P) > = 



+fK?n 2 K 



due 



-iuP-x 



m s jo 
i{Pfi^ v P u x^ 



■2P-Xj 

due~ luP ' x (p p (u) , 



B(u), 



f K Ml 
6m 



s jo 



due 



-iuP-x 



I 



f?,K {{PfxPaQvP - PvPag^p) ~ {P^P^t 

-p»pp9ia)} J va mK {aiy- iP<as+vag) 



<0\u(0)w 5 g s G a p(vx)s(x)\K(P)> = P, 



P a Xp - P X C 



P-x 



fKm 2 K 



-iP-x{a s +va g ) 



J VaiA\\(ai)t 
+fKm 2 K (P l3 g a ^ - P a 9fo) 
J Va i A ± (a i )e- iP < aa+va " ) 

P q Xr Pb x oc r 2 

P„ ^ — - — f K m K 

\ —iP-x{a s +va g ) 



P-X 



J VaiV\\{ai) 
+fKm 2 K (Ppg a ^ - P c 

Va i V ± (a i )e- iP < as+vas) 



(23) 



here the operator G a p is the dual of the G a p, G a p = \e a p llv G llv , Vcti is defined 
as Voti = daida 2 da 3 S(l — ct\ — a 2 — « 3 ), a 2 , a 3 ) = A± + V± + A\\ + V\\ and 

^(ai, a 2 , as) = A\\ + V\\ — 2A± — 2V±. The light-cone distribution amplitudes are 
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parameterized 

<f>K(u,fl) = 

<p p (u,fi) = 



as 



A ± (ai,fi) 
A(u,n) 



6«(1 - u) |l + a x c} (2m - 1) + a 2 cj (2m - 1) + a 4 c| (2m - 1)} , 
1 + { 307/3 ~\p 2 }cl (2m -l) 

+ {-3^3 - ^P 2 - ^pV} c|(2u - 1) , 



(Pa(u,fj) = 6u(l-U) \ l + 



T(ai,n) 
V\\(ai,fx) 



1 7 2 3 2 

360a u a s a^ |l + A 3 (a u - a s ) + cj 3 i(7a 9 - 3) 

120a u a s a g (v 00 + v 10 (3a g - 1)) , 

120a u a s a g a 10 (a s - ot u ) , 

-30a 2 {h 00 (l - a g ) + V - a g ) - Qa u a s 



C|(2u-1) 



10 



30a 2 (a u - a s ) j/i o + V« 9 + ^io(5a 5 - 3)| , 

, f 16 24 20 
6m(1 - m) < — + — a 2 + 20773 + yr? 4 



+ 



+ 



1 1 7 

-T5 + l6-27^ 3 



10 
27 



C|(2m-1) 



11 
210 



a 2 



4 

135 



%^3 



C|(2m-1) 



-a 2 + 21774^4 



9k(u,ii) = 
B(u,fi) = 



{2m 3 (10 - 15m + 6m 2 ) logw + 2m 3 (10 - 15m + 6m 2 ) log 
+mm(2 + 13mm)} , 

l + ^c|(2u-l)+^c|(2«-l), 
9k(u,h) - (f> K (u, fi) , 



u 



(24) 
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where 

n oo — v oo — , 
21 9 

OlO = "^4^4 - — a2 , 
o zU 

21 

o 

/, 7 3 

7 3 

18 20 
#2 = 1 + y «2 + 60r/ 3 + y 774 , 

9 

94 = -7^2-67/3^3, (25) 

- - - 4- 2 

here C 2 2 , C 4 2 and C 2 2 are Gegenbauer polynomials, 773 = j^- mq M T a and p 2 = ^f- 

[HIE1 Ell HH [13- 

The explicit expressions of the Borel transformed correlation functions Bm^X 1 
and B^II^ in the hadronic representation, 

D „! / G D * DsK m D *f D *f Ds m 2 Ds m 2 D * -m 2 Ds +m 2 K ( m 2 D , m 2 Dg \ 
BM ^ ~ I r^T^ eXP \~MT" MfJ 



G Do D 3 KfD fD s m 2 Dg r mf, m|, 



+ m c + m 5 3 6XP I" " Mf ^ + 

G D *D a Km D *f D *f Ds m 2 Ds m 2 D , + m 2 Ds - m 2 K ( m 2 D , m 



_|_ 4 : : : : ill _ m _ eX p 

1 m c + m s m 2 D * \ M\ M? 



2 



r w^.-,, fflp ._^_^np + ,„ | (26) 



'2 



G Do D aK fD fDm 2 Ds ^ { _ mj 

+ m s 6XP \ Ml M 

Km Dt f Dt f D rn 2 D m 2 Dt -rn 2 D + m 2 K J m 2 
1 m c + m M mf,. ^ I M 2 M| 



Gp*E)K m D* f D*f D m p m Dt -m D + m K m D% ml 
dmII,. = < o exp < — — ^ — 



_ G D3oDK f Ds0 f D m 2 D j _ m 2 Ds0 



m c + m u ~" r { Ml Ml 
G D * DK m D *J D *J D m 2 D m 2 Di + m 2 D -mf K j m£. m 2 , 



m c + m u m 2 -,, ' j M 2 M| 



GD s DKfD s0 fpm 2 D f m D 

m c + m u GXP \ Ml M 2 



+ " exp ^ - ^ ^ + ■ ■ ■ , (27) 




here we have not shown the contributions from the high resonances and continuum 
states explicitly for simplicity. 
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The explicit expressions of the Borel transformed correlation functions -B^n* 



and B M U. 2 at the level of quark-gluon degrees of freedom, 



q^exp 



M 2 



u (l- u )m 2 K + m 2 c \ / f K m 2 K M 2 



<f>p( U o) - ^l—^aiUo) 



puo ruo pl—a s 

+m c f K m 2 K / dtB(t) - h K m 2 K I da s 

JO JO Ju Q -Ct s 



da. 



d 
6duo 

[2(a s - u ) - a g ] T(a u , a g , a s ) 



at 



2m c f K m i K f 1 



M 2 
2m c f K m A K 



l-UQ 

1-uo 



UQ — a s 

1 - uo)$(l - a - (3, (3, a) 
a 2 



o 



M 2 

$(1 — a — a g , a g , a) 



otg pi— 13 

d/3 

Jo 

/•uq pa a pi fl-°<g f a a 

da g / da s \ da + / da g \ da s / da 

J uo—a„ Jo Jl—uo J uo—a„ Jo 



a n 



+ Py, exp 

fi<m 2 K 
6m s 

+m c f K 



«o(l - u )m 2 K + m 2 c \ / f K m 2 K M 



M 2 



1 ( f K m 2 M 2 1 . , r 2 r° , „. , 
f < u (p p {u ) +m c f K m K u / diS(i) 

J I m s Jo 



d 



b a (u )M 2 + —u (j) a (uo)M 2 + 2m;;0 (7 (w o ) 
(p K {uo)M 2 + 



m 2 K mlA(uo) 



+3u f 3K m 2 K I da 



AM 2 
da 



"0— «s 



2u m c f K m 4 K f 1 



M 2 
2u m c f K m A K 



M 2 - u m 2 K 



da ( , 



1— «0 
1— WO 



d/3 



1— a. 



da 



9 a 2 

uo-a s 
1-13 



'P\Oiui Qtgi ^s) 



da- 



1 - m )$(1 - a - f3, ft, a) 



at 



M 2 

$(1 — a — a g , a g , a) 
a„ 



pl—uo puo pct s pi pl—a g p 

I da g I da s I da + da g / da s 

JO J UQ — OLg JO J 1 — WO JUQ — Clg JO 



da 



+ m c f K m K / da 



1— a s 



da, 



ty(a u , OL g , a s ) 



Bm^I = BmII^u < — > l-u;a s < — > a u ) , 

hprP „ n - M l ., M 2 - M f M l. 



(28) 
(29) 
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AA 



The explicit expressions of the notations AA, BB, CC, DD, EE and FF, 
u (l - u )m 2 K + m; 



cxp 



M 2 



cxp 



s% \ \ S f K m 2 K M 2 



fKm 2 K 
u (l - u )m 2 K + m; 



-m c f K (f) K (u )M 2 - 
+ exp 



d 



r (u )M 2 + —u (f) a (u )M 2 + 2m 2 c (j) a (u Q ) 



M 2 

uo /•!— a s 



ruo pi—as 

+3u f 3K m 2 K / da s / da 

Jo J UQ—a s 



| ^rn c f K m 2 K u J dtB(t) 



+ 



f K m 2 K m 3 c A(u ) 



4M 2 



d 



M z - u m 2 K 

duo 



+ 2f 3 K 

2u m c f K m 4 K f 1 



M 2 
2u m c f K m 4 K 



da. 



1— no 
1— wo 



da 
d/3 



a g 

1— a» 



c!a; — r 



T(a u , a g , ck s ) 



1-/3 rfQ (!-%)$(! -a -/?,/?, a) 



at 



o 



M 2 

$(1 — a — a g , a g , a) 



da r 



uo 



da* 



u -a g 



da + 

-l—a. 



da. 



+ m. 



a. 



fi<m 2 K I da s I da g ^^ au ' a(n 

Jo Ju -a s a g 



l-OCg 

da s 

UQ—a g 

_a s ) 

9 



ra s 

I da 
Jo 

(30) 



BB 



ex P < + - exp I -M\ 



M 2 



-m c f K (f) K (u )M 
+ exp 



2 _ fKrrijc 
6m s 

u (l - u )m 2 K + m 



d 



<p a (u )M 2 + — M O ff (« o )M 2 + 2m 2 c (j) a (u ) 
du 



M 2 



| ^jn c f K m 2 K u J dtB(t) 



+ 



f K m 2 K m 3 c A(u ) 
AM 2 



ruo rl—a s 

+3u f3 K m 2 K / da s / da. 

JO J un — a<. 



T(a U) a g) a s ) 



a n 



3K 



M — u m K 



du 

2u m c f K m 4 K r 



uo-a s "-g 
uo /-l— a s 

da s 

J uo— a s 



a g + a s -u 
da,,— 



at 



T(a u ,a g ,a s ) 



M 2 
2u m c f K m 



l-u JO 



da q I d(3 



as ^ f'~\ (l-u )$(l-a-(3,(3,a) 



da 



o 



at 



M 2 



$(1 — a — a g , a g , a 



1— uo ruo ras rl fl— a, 

dct g / da s I da + J da, 

J UQ—a g Jo J 1— uo J uo—a g 

I— a. 



da s I da 



+ m. 



f K m K / da s J da g — — 

Jo Ju -a s a g 



(31) 
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M 2 



M 2 



to. 



+ exp ^(i + | mc/A . m|{ ^ dtB(t) 



-f3Km 2 K / 
2m c f K m 4 K >' ' 



M 2 
2m c f K m% 



1 ° s rfa [2(a s - u Q ) - Og\ T(a u , a g , a s ) 
9 a 2 

u -a s 

a g fl-/3 



da g / d/3 
l-uo Jo 



da 



(1 - u )$(l - a - 13,(3, a) 



at 



M 2 

$(1 — a — a g , a g , a) 



rl—uo ruo ra s rl pi— a g ra s 

/ da g I da s / da + da g / da s / da 

Jo J uq— a g Jo J 1—uq Juo—a g Jo 



a n 



(32) 



cxp 



u (l - u )m 2 K + to; 



M 2 



— exp 



'D s0 



M 2 



f K m\M 2 



TO,. 



+ exp {~ Jl™^™*./, diB(i) 



-fzKm K 



da« 



da n 



u -a s 



at 



2m c f K m 4 K f 1 
M 2 

2m c f K m% 



da n 



M 2 



l-uo 
l-uo 



d/3 



1-/3 



da 



(1 - u )$(l - a - (3, (3, a) 



at 



da, 



/ da s / 

J U{)—Ctg JO 



da + 



1— uo 



rl-a g rc 

/ da s \ 

J uo—a g Jo 



da 



$(1 — a — a g , a g , a) 



a n 



(33) 
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EE = 



u (l - u )m 2 K + m\ 
(>XP 1 M 2 <~ ~ <>xp 



XXffKrnZc^ 



-m c f K <j> K (u )M 2 - 



d 



m* 



Uo4>p(uo) 



^(u )M 2 + —u^ a {u )M 2 + 2m 2 c <t> a {uQ) 
duo 



+ exp |_ Mo (l - y + ml | { mjKmluo £> dtB(t) 



+ 



f K m 2 K m 3 c A(u ) 
AM 2 



+3u f 3K m K / da s / da g 

'0 J UQ—a s 



d 



M z - Uom 2 K 

duo 



+ 2f3K 

2u Q m c f K m 4 K f 1 



da 



a g 

1— a s 



a g + a s - u Q 

da 



s I \j,u.q „ 

a 



T(a u , a g , a s ) 



M 2 

2u m c f K m 4 K 
M 2 



da„ 

l-u 



d(3 



uo—a s 

1_/3 da^— UQ) ^ 1 - a -^^ a ) 



1-uo <>0 ./O a 



da n 



o 



«0 



da. 



$(1 — a — a g , a g , a) 



uo—a g jo 

UQ 



a n 



+ m cfK m K / da 



da + 



1 — «s 



da r 



l-uo 



l—a a 



da* 



UQ—a g 







da, 



a 



da 



(34) 



cxp 



« (1 - u )m 2 K + m 2 c 



M 2 



— exp 



d s0 1 1 f /*<M 2 



M2 J J I 



m,. 



-m c f K (p K (uo)M 2 



fi<m 2 K 
6m« 



d 



^(m )M 2 + — uqM u o)M 2 + 2m 2 (T (w o ) 
duo 



+exp p (i-„o)< + m; j | mc//fm|flIo ^ diB(i) 

r 1 ""' T(a u a , a s ) 
+3u f 3K m 2 K / da s / da 9 g ' 



+ 



f K m 2 K m*A(uo) 
AM 2 



uo—a s 



d 



M 2 - u m 2 K 

duo 



+2f 3 K 

2u m c f K m 4 K < 
M 2 

2u m c f K m 4 K 



a 

da. 



a g + a s -u 

da 



s I — g 2 



da, 



i— «0 



[ 9 d(3 [ 
Jo Jo 
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uo—a s "-g 
1-/3 

da 



M 2 

$(1 — a — a fl , a g , a) 



\ da g I da s / da + da g 

Jo Jun—a a Jo J 1 — \ 



rl-otg pc 

I da s \ 

J UQ—Clg JO 



' UQ — Oig JO Jl—UQ J UQ—Oig 

,2 



a r , 



+ m c f K m K / da 



o 
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da 
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(35) 
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